Spin-orbit coupling effect on quantum Hall ferromagnets with vanishing Zeeman 

energy 



On 
On 



O 

o 



Vladimir I. Fal'ko* and S.V. Iordanskii** 

* School of Physics and Chemistry, Lancaster University, LAI J^YB, UK 

* Landau Institute for Theoretical Physics, Kosygina 1, Moscow, Russia 

(today) 

We present the phase diagram of a ferromagnetic v = 2N + 1 quantum Hall effect liquid in 
a narrow quantum well with vanishing single-particle Zeeman splitting, ez and pronounced spin- 
orbit coupling. Upon decreasing ez, the spin-polarization field of a liquid takes, first, the easy- axis 
configuration, followed by the formation of a helical state, which affects the transport and NMR 
properties of a liquid and the form of topological defects in it. 
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The concept of quantum Hall ferromagnets (QHF)trEI 
is now broadly accepted as an approach to treat the 
2D electron systems in the vicinity of odd-integer fill- 
ing factors, v = 2M-.+ 1, in particular, due to the 
cess in the theorySo and experimental observations 
of quantum Hall effect skyrmions, which are electrically 
charged topological textures in the ferromagnetic order 
parameter|-|Ua two dimensions. The attempts to stabilize 
skyrmionsau, as compared to electrons and holes on the 
top of a filkd, spin-split Landau level (LL) have triggered 
the studiesa'B of semiconductor structures with a reduced 
value of the single-electron Zeeman energy, ez- 

The latter is possible to realize using GaAs/AlGaAs 
structures, due to a strong spin-orbit— (SO) coupling in 
this zinc-blend-type semiconductorlljliil. In particular, 
the conduction band electron Lande factor, g in GaAs 
may reduce its absolute value or even change sign under 
a hydrostatic pressure, and also due to the electron cop-, 
finement in a narrow quantum well (QW) . Both theorytHl 
and experimentlij'E-j agree on that in a non-strained 55A- 
wide GaAs/A1.33Ga.67As QW, g — 0, and the value 
g = —0.1 has rbeen recently measured in a non-strained 
w = 68A QWu. As a result, a narrow quantum well is 
a system, where Zeeman energy can be swept through 
zero, and where the tendency of interacting electrons 
at v = (2N + 1) to form a ferromagnetic state (with 
a large exchange energy 3 > Ez) polarized along the 
external -magnetic field confronts the effects of the SO 
couplingtlxliJ itself, as an alternative source for choosing 
the spin-polarization direction for 2D electrons. 

In GaAs/AlGaAs quantum wells grown along the [001]- 
facet, the effective two-dimensional SO coupling has the 
form of H so = u so {p x a x - p y o- y )+u so e l: > z piO- : > , where p = 
—iV — eA/c is the electron momentum (A = (0,a;i? z ,0)), 
<7j are Pauli matrices, and ft = 1. The first term in _ff so 
reflects the lack of inversion symmetry in the well grown 
along the [001] facet with square 2D lattice symmetry 
and comes from the (jp^u) jarm=parabolicity in the con- 
duction band of bulk GaAsiMS, sf) that u so = j(pl)- 
Experimentally^ and theoreticallyEj found values of 7 
range from 7 = 27eVA 3 to 22eVA 3 . The second SO term 



is due to the quantum well potential asymmetryfij. In a 
narrow QW, with ip — 2~ 1 / 2 su\{zit/w) transverse p^ 
of the electron wave function, u so — 7 (ir/w) ^ u s 
(Directions for x- and y-axes are chosen in such a way 
that u so ,u so > 0.) 

In the present paper, we study the effect of spin-orbit 
coupling on the v = (2 A + 1) QHF formed by electrons 
with a vanishing single-particle Zeeman energy in a nar- 
row QW in a perpendicular magnetic field. In this analy- 
sis, we address the properties of narrow wells, since they 
have a more prominent SO coupling in the 2D electron 
Hamiltonian. In a high magnetic field providing to c — 
\eB z \ jmc 3> u so /A (A = \fc/ \eB z \), a weak SO coupling 
does not alter Landau quantization. However, it affects 
the evolution of ferromagnetic properties of a quantum 
Hall effect liquid upon sweeping ez through zero. As in 
ordinary ferromagnets—. SO coupling results in the crys- 
talline anisotropy fieldta, that deflects spin polarization n 
from alignement in the magnetic field direction, \ z . As a 
result, a liquid with |ez| < e z = 2 (w so /Aaj c )~ 3jv+i takes 
one of two two- fold degenerate easy- axis magnetic states, 
which generates a source for formation of a domain struc- 
ture in a QHF. Moreover, as a feature of the SO coupling 
in the 2D electron system lacking inversion symmetry, 
for a tinier splitting, \ez\ < e| = \/2 (u so /\u c ) Sat+i, 
the spin-polarization acquires a helically twisted texture 
with a mesoscopic-scale period, £/A « 3.25uj c X/u so , in 
the [110] crystallographic direction. The transition of 
a liquid into the helical state may manifest itself in a 
change of the NMR-lineshape from the QW structure, 
or in the anisotropy of dissipative transport characteris- 
tics of a QHF. We also discuss topological defects in the 
helical state, ±e/2-charged dislocations, which, in pairs, 
constitute skyrmions. 

The recipedu of describing smoothly textured QHF's 
at odd-integer filling factors is to use the 2D sigma-model, 
which operates with the energy functional, $ {n(x)} of 
2D electrons expressed in terms of their local excess spin 
polarization field n(r,t) (|n| =1). Locally, the QHF can 
be viewed as a liquid of electrons which fully occupy 
(N + 1) LL's with spin parallel to n, and N LL's with 



antiparallel spins. This assumes the existence of a lo- 
cal unitary spin-transformation, U(r) of electron wave 
functions, which reduces the filling of the (N + l)-st LL 
to a complete occupation of only states 'up' with re- 
spect to the locally determined axis n(r), and which is 
related to the spin-polarization field as n % — Tr(s z A), 
where s'(r) = U{v)a l U\v), and A = {l + a z )/2 is 
the electron spin-density matrix in the rotated frame. 
Such a liquid retains incompressibility, which is guaran- 
teed by a large exchange energy gap, 3. The deriva- 
tion of a sigma-model consists in the use of Hubbard- 
Stratonovich transformation and the saddle-point self- 
consistency equation, as a method to obtain an expan- 
sion of the thermodynamic potential of a liquid over small 
gradients of a polarization field, or, equivalently, over 
the matrix Tt(r) = U(r) (-iW^r)) = (iW(r)) E/t(r). 
The latter matrix appears in a local perturbation to the 
single-particle Hamiltonian written in the rotated spin- 
frame, 
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where a T = ±(ip x ^f ap y ) / \^2 are the inter-LL operators, 
a = eB z / \eB z \ indicates the direction of the cyclotron 
rotation of carriers, and fl± = (Cl x =F iafl y )/v2- Further 
procedure (equivalent to the Hartree-Fock approxima- 
tion) consists in a perturbative expansion of the saddle- 
point Hubbard-Stratonovich action (obtained for a given 
single-particle Hamiltonian, which includes H£-t + H so ) up 
to the second order both in £1 and the SO coupling, H so 
(assuming that u so / A < 3). In the local spin- frame, 

H so = (ugo/X) [i(a+s + - s~d-) + n + s + + s~fi_] 

+ (au so /A) [a+s - + s + a_ — i(fl + s~ — s + J7_)] , 

where s ± = (s x ^fais y )/y/2, and the relevant part of H so 
is off-diagonal with respect to the LL-numberc3. 

The calculation, which leads us to the sheet density 
of thermodynamic potential of a QHEJ>{n(x)} differs 
from earlier Hartree-Fock calculationsQO only by taking 
into account the SO-coupling term H so , along with Hq. 
It results in 
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where nil is the planar component of the spin-polarization field, and n x — n y , n v — n x . To obtain <I> in Eq. p), we 
have extended perturbative expansion up to the terms 3V 2 , 9 (u so /Xuj c ) V and 9 (u so /\u> c ) . The exchange-factor, 
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was calculated for each odd-integer filling v = 2./V + 1 (N = 0, 1, 2, ...); L^(z) are the generalized Laguerre polynomials. 
For V(r) = e 2 /r\, which is a reasonable approximation for the 2D electron interaction in a narrow QW (w <C A), 
6\ = i, 8-2 = jg, 63 ~ 0. 57, 84 « 0. 67, and #5 « 0. 76. Since the relevant part of the SO coupling is off-diagonal with 
respect to the Landau level number, it appears only in the second-order of a perturbation theory, or due to its mixing 
with the iJo-termEl. The single-particle spin-splitting, ez in Eq. (|I]) is corrected by the effect of the 2D single-particle 
SO-coupling: ez — pgB — av (u 2 — w 2 ) /w c A 2 . We also include into $(n, Vn) the topological term and the Coulomb 
energy of additional charges, p(x), in order to discuss charged skyrmion-type textures, 
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From a phenomenological point of view, thermody- 
namic potential $ {n(x)} in Eq. (fy) describes an easy- 
axis ferromagnet with square 2D Bravais lattice and bro- 
ken inversion symmetry, and in a perpendicular magnetic 
field. It contains all terms in the magnetization energy 
allowed by the crystalline symmetry of the [OOlJ-grown 
quantum well in a zinc-blend- type semiconductorE3. The 
sketched above microscopic derivation of ${n(x)} was 
necessary to obtain the values of coefficients in front of 
the phenomenologically allowed invariants composed of n 
and V. The first term in $ describes spin-stiffness. The 
second term determines an easy-axis anisotropy along 
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1+ = [110] J \[2. For structures with m so ^> u so , it rather 
defines a weakly anisotropic easy plane for spin polar- 
ization, which competes with the Zeeman energy term. 
Such a competition resumes itself in the deviation of po- 
larization from a fully ^-aligned state at \ez\ < £7 
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As a function of a varying ez, this can be viewed as a 
second order phase transition into the easy-axis state, 
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across which the symmetry between ±1+ magnetic direc- 
tions gets spontaneously brokenEj. 

Since for \ez\ < s\ both easy-axis configurations, n + 
and n_ have the same energy density, 



*o 



%N+1 (u+/\ 

'4vrA 2 



1+(S 



(5) 



the easy-axis state of a QHF tends to acquire a do- 
main structure: by splitting dynamically into the set of 
mesoscopic-size regions with opposite polarization pro- 
jections onto the [110]-axis. The latter possibility has 
to affect the skyrmion-dominated dissipative transport 
properties of a liquid. The matter is that the activation 
energy of a skyrmion-antiskyrmion pair confined to the 
domain wall is lower, than in the 2D bulkcil. In fact, the 
larger is the difference in the polarization between two 
domains, the more skyrmion is energetically confined to 
it. Since n±-polarized states are degenerate, the areas 
covered by n + and n_ domains are statistically equal, so 
that the network of better conducting domain walls (with 
a lower activation energy for thermally excited carriers) 
forms a percolation cluster, thus resulting in a continu- 
ous decline in the activation energy for the macroscopic 
a xx which would follow the decrease of £z- 

A further analysis of the functional in Eq. (fy) extended 
onto the limit of Ez = shows that, apart from the do- 
main structure formation, there is another reason fcp-jthe 
field nfr) to be inhomogeneous. The fourth termca in 
Eq. (|l|) tends to twist the polarization field of a QHF 
into helical texture, 



n(r) = lsin(0(r)) + 1 2 cos(^(r)). 



(6) 



The latter is characterized by the helicity plane built 
upon two unit vectors, 1 = (l x ,l y ,0) and 1 2 , spatial ori- 
entation m, and period £, <j)(x + m£) —<f>(r) + 2tt. As a 
variational approximation, one can use n(r) in Eq. (J6|) 
with </>(r) = mr/£, treating 1, m and C as minimization 
parameters. The resulting texture can be viewed as an 
image of a spoke in a wheel rolling in the direction of 
m, with <f> being the integral angle encircled by a spoke, 
so that we shall call <f> the 'helicity phase'. The energy 
density of an optimal variational state, 
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is lower than the energy Q c (ez — 0) of a homogenous 
easy- axis configuration in Eq. (H). Optimal variational 
state provides us with values of $h and C which are vepy, 
close to the n-field distribution that really minimizesE3 
the functional $ in Eq. (|l|). For the experimentally 
relevant example of u + ~ u_ (i.e., u so ^> u so ), <I>^ = 
-§ (wso/Awc) 2 (3jv+i/47tA 2 ) and C v in Eq. (@) should be 
comparedE3 with $ h « -1.53 (u so /Xlu c ) 2 (^ n+ i/4wX 2 ) 



and C « 3.25w c A 2 /m so . We, therefore, assess the stabil- 
ity of a helical state, against the easy-axis one, on the 
basis of the energetics of variational helical texture with 
parameters given by Eq. (|7|). 

Due to the difference in the symmetry of a helical 
and easy-axis states, which cannot be continuously trans- 
formed one into another, the transformation between 
them can be viewed as a first-order phase transition. The 
condition for such a transition, <&h = < l ) c(£z), determines 
critical value of Zeeman splitting, e\ estimated as 
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The helical state formation can manifest itself in sev- 
eral observations. For example, the local value of the 
Knight shiftH, 5\^ in the spin-splitting of Ga and As nu- 
clei located in the QW acquires an alternate coordinate- 
dependent sign, thus modifying the NMR lineshape, 1(5). 
Locally, the NMR shift 6 = (ui — iv ) /5hi is due to the 
hyperfine interaction of nuclear spins with fully polarized 
electrons, with Shi being its maximal value just in the 
QW center. In a homogeneously polarized gas, tie NMR 
line from the QW has a double-peak structurecl, 1(1 > 

I/O 

5 > 0) oc [5(1 — 5)] , with a distinct satellite at 6 = 1 
split by the hyperfine coupling. In the easy-axis con- 
figuration, Eq. (H), the double-peak structure persists, 
with a reduced maximal splitting: <5hf — > (sz/sz,) 5m, 
as far as n z = — £z/s|. On the contrary, in a helical 
phase, this has to transform into a single broadly tailed 
resonance with a non-Lorenzian shape approximated by 
1(5) oc \5\' 1/2 for 1 > 5 > -1. 

The anisotropy of transport characteristics of a QHF 
with respect to [110] and [110] crystallographic direc- 
tions may be another feature of the helical state. Speak- 
ing about dissipative conductivity formed by thermally 
activated electron- hole pairs at the spin-split LL's, this 
can be understood after taking into account that charge- 
carrying excitations determined in a locally rotated spin- 
frame (adjusted to n(r)) are subjected to a smooth 
(C ^> A) potential due to Zeeman energy, (ez/2) cos([x — 
y\/y2C). Therefore, at low temperatures, T < ez the 
dissipative conductivity er_._ along m ||[110] would be 
suppressed, as compared to o+ !+ (across m). 

For the dissipative transport dominated by 
skyrmionstj, the difference between u~- and u+,+ is 
to be the result of the anisotropy of a skyrmion itself. In 
fact, the form of a skyrmion in a helically twisted texture 
is quite complex: In a periodic system, these are dislo- 
cations which represent the very elementary topological 
defects, rather than skyrmions. The periodicity of a he- 
lical texture in Eq. (6) is controlled by the helicity phase 
</>o(mr/£) in Eqs. (I) and (M). One missing (or extra) 
period in one half of a plane, as compared to the other 
half (a dislocation) is equivalent to the phase shift of 
±2-7r accumulated at large distances from the dislocation 
core, thus resulting in the winding number D = ±1. On 
the other hand, we assume that the dislocation core is 



not singular. To illustrate the topology of a non-singular 
core, let us draw a large-radius circle around a disloca- 
tion. At large distances, where helical structure is not 
perturbed, such a contour maps into the equator of a 
unit sphere and encircles it N or (N — 1) times, depend- 
ing on which one of two semi-circles is retraced: drawn 
above, or below the dislocation. 

Upon moving the upper half of a contour down through 
the dislocation, an extra loop encircling the unit sphere 
equator should continuously disappear. The latter is 
possible if the contour image slips through either the 
+1+, or, alternatively, — 1 + pole of a unit sphere, which 
can be modelled by such a field configuration n(r) = 
\ z n z + l + n + + l_n _ , that 
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n + (r — > oo) = 0, 



where (r, ip) are polar coordinates calculated from the dis- 
location center, r = 0. The image of a 2D plane provided 
by n(r) maps onto only one half of a sphere |n| = 1, so 
that the dislocation core is characterized by an additional 
topological number, $ = n + (0) = ±1 distinguishing be- 
tween 'left' and 'right' semi-spheres. Using /?(x) in Eq. 
(H), we find that the core of a dislocation (D = 1) or anti- 
dislocation (D = —1) carry a half-integer electric charge 
Jofxp(x) = \"3D. However, an isolated dislocation has a 
logarithmically large energy, 
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and, at low temperatures, dislocations and anti- 
dislocations have to form pairs bound by a long-range 
(logarithmic) attraction, except, maybe, for a possible 
Kosterlitz-Thouless melting effect. Since both disloca- 
tion and anti-dislocation in a bound pair may be equally 
charged, such a pair, (D, $) and (-D, — i?) together con- 
stitute a skyrmion. 

The result of the above analysis of phases of a ferro- 
magnetic quantum Hall effect liquid in a narrow QW can 
be summarized as follows. Upon decreasing the single- 
particle Zeeman splitting, \ez\ (e.g., by pressure), spin 
polarization of a liquid starts to acquire at \sz\ = £z the 
easy-axis configuration, which is followed by the abrupt 
fall into a helical state at \ez\ = £% « e|/V2. Us- 
ing the bulk SO coupling parameter 7 = 25eVA 3 , as 
a reference, we estimate for the v = 1 liquid in a 68A- 
wide GaAs/AlGaAs quantum well structure with .a car- 
rier density 2.8 x 10 11 cm -2 studied by Maude et a!3 that 
e z /(e 2 /xA) « 1.5 x 1CT 3 , and e\/{e 2 / X \) ~ 1 X 1(T 3 , 
which roughlypfits into the range of a variable Zeeman 
energy in Ref.fi, where the dissipative transport activa- 
tion has been drastically affected by pressure. The heli- 
cal texture period estimated for the same parameters is 
£ - 5 x 10 3 A. Note that, according to Eqs. (|) and (|), 
the parametric range of pressures where helical and easy- 
axis phases are stable is broader for higher odd-integer 
filling factors in the same density structure. 
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